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XXIV. — On the development of the disturbing Function, upon which depend 
the inequalities of the motions of the Planets, caused by their mutual attrac- 
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THE perturbations of the planets caused by tbeit annual at.raet.on depend 
chiefly upon one algebraic expression, from the development of which all the 
inequalities of their motions are derived. This function is very complicated, 
and requires much labour and many tedious operations to expand it in a series 
of parts which can be separately computed according to the occasions of the 
astronomer. The progress of physical astronomy has undoubtedly been re- 
tarded by the excessive length and irksomeness attending the arithmetical 
calculation of the inequalities. On this subject astronomers generally and 
continually complain ; and that their complaints are well founded, is very aptly 
illustrated by a paper contained in the last year's Transactions of this Society. 
The disturbing function is usually expanded in parts arranged according to 
the powers and products of the excentricities and the inclinations of the orbits 
to the ecliptic ; and, as these elements are always small, the resulting series 
decreases in every case with great rapidity. No difficulty would therefore be 
found in this research, if an inequality depended solely on the quantity of the 
coefficient of its argument in the expanded function ; because the terms of the 
series decrease so fast, that all of them, except those of the first order, or^ at 
mostj those of the first and second orders, might be safely neglected, as pro- 
ducing no sensible variation in the planet's motion. But the magnitude of an 
inequality depends upon the length of its period, as well as upon the coefficient 
of its argument. When the former embraces a course of many years, the latter, 
although almost evanescent in the differential equation, acquires a great mul- 
tiplier in the process of integration, and thus comes to have a sensible effect 
on the place of the planet. Such is the origin of some of the most remarkable 
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of the planetary irregularities, and in particular, of the great equations in the 
mean motions of Jupiter and Saturn, the discovery of which does so much 
honour to the sagacity of Laplace. It is not, therefore, enough to calculate 
the terms of the first order, or of the first and second orders, in the expansion 
of the disturbing function. This is already done in most of the books that 
treat of physical astronomy with all the care and fulness which the importance 
of the subject demands, leaving little room for further improvement. In the 
present state of the theory of the planetary motions, it is requisite that the 
astronomer have it in his power to compute any term in the expansion of the 
disturbing function below the sixth order ; since it has been found that there 
are inequalities depending upon terms of the fifth order, which have a sensible 
effect on the motions of some of the planets. 

A research that has for its object the lessening of the difficulty attending 
the expansion of the disturbing function, and the bringing of that expression 
more under the power of the astronomer, is one of considerable interest, and of 
some moment to the progress of physical astronomy. But the question is not 
to exhibit separately every individual argument with its coefficient, which 
would be of little utility in practice, since, although their number be infinite, 
very few of them are of any account in computing the place of a planet. The 
end to be aimed at, is to give the function such a form that the astronomer 
may have it in his power to select any inequality he may wish to examine, and 
to compute the coefficient of its argument by an arithmetical process of mo- 
derate length. The present Paper is an attempt of this kind. The investiga- 
tion comprehends every argument not passing the fifth order ; but, as the for- 
mulas are regular, the method may be extended indefinitely to any order. 

1. Let x, y, z, represent the rectangular coordinates of the disturbed planet, 
z being perpendicular to a fixed plane passing through the sun's centre, 
at which point the origin of the coordinates is placed: in like manner let 
x', y, z f be the coordinates of the disturbing planet ; and put m f for the pro- 
portion of the mass of this planet to the sum of the masses of the sun and the 
disturbed planet ; then, R denoting the disturbing function, we shall have, 

If r and r f denote the distances of the disturbed and disturbing planets from 
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the sun's centre ; *>, their angular distance seen at the same point ; and g, their 
rectilinear distance from one another, then, 

g> = fjr 2 — 2rr ! cos co + r 12 
R 1 r cos a) 

It appears, therefore, that the function to be expanded consists of two parts, 

of which — , the reciprocal of the distance of the two planets, is common to the 

disturbing functions of both ; but this is not the case with the other part, in 
which r and r' do not enter alike. For this reason it will contribute to distinct- 
ness to expand the two parts separately. But, previously, it is necessary to express 
cos co in terms of the angular motions of the planets in their respective orbits. 

If v represent the angular distance of the disturbed planet in its orbit from 
a fixed origin, and P the place of the node, that is, of the intersection of the 
orbit with the immoveable plane of xy, P being reckoned in the same plane 
and from the same origin as v 9 then v — P will be the angular distance of the 
planet from the node. Further, if A be the celestial arc between the node and 
the intersection of the orbits of the two planets, the distance of the planet 
from the same intersection will be equal to 

v - P - A. 

If the foregoing symbols be accented and transferred with like significations 
to the disturbing planet, the celestial arc between that planet and the inter- 
section of the orbits will be equal to 

v'-F- A f . 

Now v — P — A and d — P f — A! are two sides of a triangle of which the arc 
m is the third side ; wherefore if I be the inclination of the two orbits, we shall 

have 

cos &) = cos (u — P — A) cos (V — P' — A') 

+ cos I sin (v — P — A) sin (t/ — - P' — A'), 
or, which is the same thing, 

cos m = cos 2 -q I . cos (y — v* — - P + F — A + A ; ) 
+ sin 2 — I . cos (y + v — P — P ; — A — A f ). 
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Let N be the place of the node of the disturbed planet reckoned in the im- 
moveable plane of x y, and i the inclination of the orbit to the same plane : 
then, if v represent the longitude of the intersection of the orbits of the two 
planets, A will be the hypothenuse of a right-angled triangle, of which i is one 
angle, and v — N the side adjacent to i % wherefore 

tan A ~ - ^j ; 
and, by the usual methods, 

A = v — N •+• tan 2 -^ sin 2 (j> — N) + -^ tan 2 y sin 4 (y — N), &c. ; 
and if we put 

* « » 

| = tan 2 -g* sin 2 (v — N) + -^ tan 2 ^- sin 4 (V — N), &c, 

we shall have 

P -{-• A = y + P — N + f« 

Let [P — N] be the value of P — N at some given epoch : then the differen- 
tials of N and P being rfN and cos i . ^N, we shall have generally, 

P - N = [P - N] -/rfN +ycos*\rfN = [P - N] - 2/kn^N, 

the integral being taken for the time elapsed from the epoch. It thus appears 
that P — Nmay be considered as an invariable arc, or as one subject to an 
almost insensible secular equation. In order to abridge expressions, let us now 
put 

« = [P - N] - 2/Wi rfN + 5 ; 

vlJiv/JU. 

^ — p — A=ZV — M — V. 

If the symbols which stand for the elements of the disturbed planet be as- 
sumed, when accented, to represent the like elements of the disturbing planet, 
we shall have similarly, 

a 1 = [F - NQ - 2/kn 2 -J rfN 1 + r, 
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The value of cos a will now be thus expressed, 

cos of = cos 2 -£- 1 . cos (v — 1/ — a + a') 

1 
+ sin 2 tt I . cos (y + 1/ — 05 — a' — 2 v). 

It remains to determine I, the inclination of the orbits, and v, the longitude 
of their intersection. Now I is the vertical angle of a triangle of which N — N' 
is the base, and i and i ! the adjacent angles, one being interior and the other 
exterior to the triangle : wherefore 

cos I = cos i cos i 1 + sin i sin i 1 cos (N — N'). 

In the same triangle, a perpendicular being let fall upon the base from the 
vertical angle, v — N and v — - N' are the arcs between the perpendicular and 
the extremities of the base ; and hence, 

/ N + N*\ sin (i - i 1 ) N - N' 

COt I V ~ / = ~^~-7^—~^R . COt 



sin(i + fO 2 



In the method here followed, the expression of cos & is as simple and as little 
troublesome in calculation, as it would be if one of the orbits were adopted 
for the immoveable plane. 

2. Let £, s, e, zs represent the mean motion, the epoch, the excentricity, and 
the place of the perihelion, of the disturbed planet at the time for which we 
are computing : if this time be near the epoch for which the elements have 
been assigned, the values at the epoch may be taken as the true values ; and 
if a great interval has elapsed, the values at the epoch may be corrected by 
their secular equations. Further, put \h for the mean anomaly, g for the 
equation of the centre, a for the mean distance, and a (1— s) for the radius 
vector : then, 

7 = a (1 - s) 9 

And the values of a and s, found by the solution of Kepler's problem, are 
expressed by the usual series, viz. 
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= e . 2 sin p 



f 5 sin 2 ju, 

,6 



"T" e 



. ,, 13 sin 3 ju. — 8 sin//, 

+ 6* . 



= e cos (h 

. o cos 2 a — 1 
+ <? 2 . <F 



2 



12 



+ e 4 . 



103 sin 4 ju, ~- 44 sin 2 p 
96 



, •„ 3 cos 3 a — 3 cos a 
4 cos 4 p — cos 2 ju, 



5 1Q97 s ^ n £jjz ^^ sni £fjjlj2 s * n ^ j« 5 125cos5ju,~ 13ocos3ju,+ lOcos^ 



960 



384 



Using always the same symbols which stand for the elements of the disturbed 
planet, when accented, to represent the like elements of the disturbing planet, 
and introducing the new characters <p and <p r in order to shorten expressions, 
we shall have, 

t/ _ a' = £' + s' - a '- + a* = p' + </. 

If these values, as well the values of r and r f , be substituted in that part of the 
disturbing function which is more easily dealt with, we shall obtain, 

= JS • | COS 2 -g- 1 . COS — <p ! + (J — </) . Tj^-jxa 



r cos a) 



r 



,/2 



«7 



or 



r cos u) 



+ sin 2 -g- 1 . cos (?) + p' + <r + a* - 2 *) . p~ryj« j 

a 1 9 * T /^ /\ (1 — 5)C0S(<T — 0-') 

Z2 . j cos 2 y I . cos (<p - ?/) . * — ^jz~y — 
- cos 2 ^ I . sin (?) - fl'). ^j^y^ 

. . o 1 T / , / ~ \ ( 1 — ^) COS (<r -f </) 

+ sin 2 T I . cos (9 + <p' - 2 v) . i ^_^ .JT_J 



1 



sin 2 -5- 1 . cos (<p + <p f — 2 v) . 



( 1 — s) sin (0- -f <r') 



(1-0 



2 



The inspecting of this formula is sufficient to show that the expansion of the 
part multiplied by sin 2 ■§ I is deducible without calculation from the expansion 
of the other part. For, as <r' is a series of the sines, and / a series of the cosines, 
of the multiples of p, the former will change its sign with jm\ while the latter 
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will not vary : wherefore the expansion of the part multiplied by sin 2 § I will 
be obtained from the expansion of the other part, merely by making (/ nega- 
tive, and writing <p + <p' — 2 v in place of <p — <p' in all the arguments. 

We have therefore only to expand the two first lines of the expression, for 
which purpose they may be thus written : 

/ /\ f /* n coscr' . (sin <r f ) 1 

cos (<p — <p ) . I (1 — *) cos <r X (HTpp + C 1 — *) sm * + ( \ Z /y j 

(~ COS 0"' f sill (T^ *) 

+ sin (p - <p') . j - (1 - *) sin* X (T^T)"* + C 1 — *) cos * X (i~=^?)» J 

Neglecting quantities above the fifth order with respect to the excentricities, 
we have, 

(1 — s) cos a = 1 — s — 



X 2 5 cr 2 (T 4 

£ + g + 24 — 


5<r 4 

£4 


cr 3 . s <r 3 ■ <r 5 





(1 — $)sin<r =* — s* — ~q+~q~ + 120 
And, if we assume, 

\ A (0) , aW » A^ 9 « i a< 5 > r 

(1 — ^) cos<r = A +A ^cos^ + A tf 2 cos2^. . . + A £°cos5^, 

(1 — s) sin <r = B (1) e sin/* + B C <? 2 sin 2 ^ B ( e 5 sin 5 \h\ 

we shall readily find, 



A = 


l 





f 

% 





e 4 
64 






A W = 




• 1 





3 

8 


e 2 + 


5 
96 


e 4 



,(0 ., 3_ 5 

8 e "+"96 



B W = 2-|-e 2 + ^ 



n (2) — J- 1 2 

A — 2 3 B "24 96 6 

„( 4 ) 29 
B =64 

B (s) - 77 



a< 3) 

A = 


3 
8 


13 

32 


e 2 


A <4) = 


67 
129 






A < ! >il 

A 48 
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Andy by proceeding similarly, 



coscr 



c i - tf) 



l\2 



( ) I (1) / t I ( 2 ) 'tO « f I (5) IK 

a + a . e cos (/J + a . e 2 cos 2 f«/ . . . + a . £ 5 cos 5 



,(2) 



(5) 



fT^ZTF = * ^ s * n V* ^ ^ sin 2 ft/ • ... + h . e f5 sin 5 f& 



4 (o) 



1 



e 



n e f4 



% 



6 (1) 2-|e« + ^ 



a 



(i) 



^ 



f O 



p'2 .4- — /j 1 



/4 



£4 






41 



12 



£4 



/2 



« 



a 



(2) _ _7_ 10 

* £ """ 3 

- li ^5 

- ^ _ _ 

1603 



e 



to 



(3)31 161 



(3) 



/2 



& 



6 ^ £4 

(4) _ 761 

— 96 



re 



to 



a 



(4 



192 



h (5) _ 2827 
~ £40* 



a 



(5) _ J49 
— 1£ 



We can now assign any term of the expansion we are considering. Let 
and i be any two positive numbers, zero included, of which the sum does no 
exceed 5 : the part of the expansion multiplied by e i e w cos 2 \ I will be, 



cos 



(<p — p') . < A a cos ^ ^ cos / f-c/ + B o sin « p sin i' \>J S- 



«A^„,.. ( 



+ sin (jp — - <p f ) . < — ■ B & sin i \h cos i ; fe/ + b A sin l ; ^ cos i (& > % 



and, by reduction, 



A a + B 6 B a -f 6 A 



• cos (<p — <p f — i p + ^ fO 



A a -f B ' • # & H- 



J„„ 



4 



A a — B b 15 a 



- — - I . cos (p 






4 



- b K J AT\ 

" 4 I . CO 



s (p — ^ f — i^ — i'j&e/) 



„ika — ±50 ,15a — o A t , / , * s .# #* 

| ^ -_. > + ^^—^^^^ J , COS (if — ^ + I f& + ?>% 
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The other part of the expansion, that multiplied by eW e'W sin 2 ^1, will be dif- 
ferent in no respect from the part just computed except in the arguments of 
the cosines, which must be changed according to the rule already laid down ; 
so that, instead of the four cosines above, these which follow must be substituted 
respectively, 

cos (p + p f — 2 v — i ft — i ft) 

cos (p + p' ' — 2 v + i ft + i l ft 1 ) 
cos (p + P' "-- 2 y ~ i ft* + {' ^') 

cos (p + <p' — - 2 v + ift — i ft). 

Every part of the expansion being comprehended in the formula, it follows 
that all the arguments will contain the mean motions of the planets, excepting 
the particular case when i = 1, i ! = 1 : for, in this case, p — ft and p f — ^'are 
independent of the mean motions; and the first cosines of the two parts multi- 
plied by e e' will be, 

COS (w ~ 7$ ! — CC -j- ®<) 

COS (w + zf — 2 y — a — - cd) . 
These cosines have the same coefficient, viz. 

A (I) « (1) +B (1) Z, (1) B (1) a 0) + A (1) i (I) 

4 4 

And this will be found equal to zero, when the values of the symbols are sub- 
stituted # . It thus appears, as far as the calculation has been carried, that 

<y> f*f)Ci (t) 

there are no terms in the expansion of —7*— except such as contain the mean 

motions and are periodical. 

1 
3. We next proceed to the expansion of -— , which is the more difficult part 

of the problem. The value of cos &, already investigated, may be thus written, 

cos a/ = cos (<p — p ! + <r — d) 

— sin 2 § I . < cos (p — p f + <r — tf) — cos (p + P* + ff + ^ ~ 2 y) 1 ; 

* The coefficient is equal to zero, because cr' = Z> : and for the same reason the coefficient 

A^ fl W - B (1) 5 (1) B^ l) . (1) -~^ (1) A (1) 
4 + 

is equal to zero. 

4 D 2 
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wherefore, 

D 2 = r 2 + r f2 — 2 r r f cos (<p — <p f + * — *% 

<£ = COS (<p — 9 + fir — <r') — COS (<p + ¥ + * + ^ — 2 *), 
f 2 = r 2 -f r f2 ~2r r cos a; = D 2 + 2 sin 2 § I . r r 1 q s 

consequently, 

— = p— sin 2 1 1. -j^g +4sin4|I.-^~ ? 2 ; 

and the three parts of — must be expanded separately. 

Expansion o/'-jy. 
By substituting a (1 — s) and d (1 — s f ) for r and r f , we get 

-L 1 

D = ^/a 2 (i - s) 2 + a' 2 (1 - s') s -2aa'(l-5) (I -. *') cos~(^^ +^ - <r') ; 

and it will be found that the following equation, in partial differentials, is true,. 



viz, 



,1 i 



ID D 

-35 = -*r C 1 - *) + -jrC 1 -*')- 

All the possible values of -jy in this equation are comprehended in the for- 
mula, 

^ = G X F; 

provided F and G verify these equations, 

G=ll(l-*) + ^(l— 

O=77(i«*)+|7(l--0. 

We may adopt for G any particular expression which verifies its proper 
equation ; and then F will be determined by the nature of the quantity sought. 
Any function of (1 — s) and (1 — /) of — 1 dimension will verify the equation 

for G ; and, taking the most simple case, G = -. , we shall have 

1 _ I „ 
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Now let 

Do = *J a 2 + d 2 — 2 aa! cos(<p-- <p' + <r — o^)* 

then first, if we suppose s = 0,s' = 0, we get 

1 l 



D -D 



and again, if we suppose s = s f , the result will be, 



l 



and thus we learn that 



D 1 — - 5 D 



F--L 



both when s = 0, s* = 0, and when $ == s'. From its equation we know that F 
must be a function of no dimensions with respect to (1 — s) and (1 — /); and 
this condition, as well as what has been shown to take place on the two fore- 
going suppositions, will be fulfilled by making F a function of Yl — f£~) 

s t _ s i 

= | ~ - : wherefore the expression of -jj will be as follows : 

] 1 /J — s\. 

— — ■ iii I. j *f- I ■- ■•■ — ■ — — — 1 

/ being the mark of a function which, in the present case, must be equal to 
jsp when -= = 0. 

Uq 1 s 

Since $ and s' are small quantities, the expression of jy may be expanded. 
By neglecting quantities above the fifth order with respect to s and s\ 

(0) 

_L — J*— _L R (1) g^ZLJ i T*< 2) (** - *)* R ( 5 ) (^ - *) 5 

D ~~" 1 -s" 1 " 13 , (l-s) 3 " t " 13 *(i-s) 3 n •(i- s )6- 

In this form it is obvious that all the coefficients may be readily found. For, 
by differentiating successively with respect to s', and supposing s = 0, $' = 0, 
in all the results, we obtain, 

•r(°) _ _L n W _ 1? n& - l dd ® * 
15 — D > ^ - d/> ** — T&\ ds /r > &c - 

But the expression is complicated ; and in order to reduce it to a more ma- 
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nageable form, every separate term may be resolved into a series of the powers 
and products of s and s f : which being done, and the quantities of the same 
order classed together, we shall find, 

1 (0) (1) (2) (3) 

~=B +B . s' +B V +B V 3 +&c. ; 

B (1) -B (0) ).,-(B (2) --B (1) ).2„' _(b (S) -B (2) ).3^ a - 



+ (b (2) -2B (I) + B (0) )/+(b (3) -2B (2) +B (1) ).3 



rfs* 



s* 



/ (3) (2) (1) (0)\ 

~(B ~3B +3B ~-B l ). 
Or, in the usual notation of finite differences, 

I (0) (1) (2) (3) (4) 

^=BT +B \s ! +B J* +B K V s +B .V 4 +&c.; 

(0) (1) . (2) (S) 

-AB .*-AB .2^-AB .3^ -AB .4s f3 s 
+ A 2 B .*> +A 2 B l \3*V +A 2 B .6*V 

~A 3 B°.* 3 -A 3 B (1) .4^6' 3 

(°) 
+ A*B \s\ 

i 
The value of jy seems now to be reduced to its most simple form as far as s 

and s ! are concerned. If it be observed that, in every term, the exponent of 
the power of s f is the same with the numeral affix of B, and the exponent of ^ 
the same with the index of A, it will appear that the whole expansion is com- 
prehended in the formula, 

A B . (is' s, 

•/ • •# i • 

i i % -f" i 

|3 being the coefficient of s l s in (rf — s) , and z', i 9 representing the several 
positive numbers, zero included, of which the sum (i ! + i) does not exceed 5, 
when the expansion is limited to quantities of the fifth order. 

(°) (0 

4. The values of the quantities B , B , &c, are next to be investigated. 
If, for the sake of abridging, we put, 

Q — a n ( 1 — /) — a a 1 (1 — s) cos (<p — <p' + <r — ?% 
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we shall get by repeated differentiations., 





d h 


Q 




d s f 


~ D 3 > 
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I 
d 2 — 
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1 dQ 1 


1.2 * 


ds 12 " 


~ 2 ' ds r • D 3 


1 


1 

d 3 -— 
a D 


3 rfQ Q 


1.2.3 


• ds 13 " 


- 2 • ds' * D 5 


1 


1 
rf D 


3 /dQ\* 
" 8 • Vds'/ # 


1.2.3.4 ' 


1 


1 

d 5 -~- 
D 


_i5(dQ\* Q 



j_2 9! 

"r" 2 • D 5 ' 



, 5 <y 



15 dQ Q2 35 Q^ 
4 * ds • D 7 +¥ ' I) 9 ' 



35 dQ Q 3 63 <y 

172.3.4.5 * tfs' 5 ~8Ws/ ' D 7 + 4 " dj • IF + 8 * D"' 

(1) (2) 

Now these expressions become equal to B , B , &c., when s f = 0^ = 0^ that 
is 5 when 

•^ ro ? / r i ?\ Do 2 + [a! % — « 2 ) 

Q — d 2 - a a' cos (<p — <p r + <r — <r ! ) = --^--^-V - "— ^^ 



Let us now put 

then, 

n_ p o s +/ 19 - /-v. L±*. 

and, by substituting these values in the foregoing expressions, and reducing^ 
we get 

„(°) i 
B = 
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Z-h f 



4 * D 3 ' 



P (2) _ 2 / J. . £\ , 2 

a ~ 8 VD ■+" D 6 / "+" 

R (3) - 5 f 1 J ./'U 9 - 6i ^J.^ , i 

*> - 16 VD + D 7 + 16 VD 3 + WJ 

(4) _ 35 /I /« \ 20- 15 A/ / /^\ , 6A 3 - 48 h + 51 / 
15 — 128 \D + D 7 "* 32™ VD 8 + D 7 + 

(5) _ 63 / 1 ,J^_\ , 175-140A // , / 4 \ , 30^-150^+135 // 2 / 3 \ 
B — ■ 256 \D + D '7 H 256 VD 3 + D V H 188 ' W + A>7 ' 



64 • D 5 



,1 I 

it — - /?// 



These quantities have been deduced from the differentials -~jj, y^ • "J" 71 



5 



1 1 

d B 1 d ^D 



&c. If the differentials -^- , -j— ^ . -yr, &c, be used, the same results will 

be found, with this difference, that f and h, retaining the same numerical 
values, will both change their signs. 

If we write «*// for (p — <p' + $ — # f ) and differentiate the formula, 



c*- 1 



D n _ & a f sin $ 



the result will be, 

dd~ 
D , a a! cos $ , , a s g^ sin 3 ^ 

now, 

a d cos ^ = —a~ — -^- , 



a f2 « 2 sinH = - ^T +^4^ • D 2 - ^— 5LL 



wherefore, 
i 
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whence we deduce 



i 

dd — 



P _ 4 ^ D n _ n 

-pn-l-4 ~~" rc . ?* + 2 * dty 2 w + 2 



J -'o ■ L 'o 



In this formula put n successively equal to 1, 3, 5, 7; then 

D 5 ~~ 3 • rf^ 2 3 * D + 3 * I) 3 ' 

/* _ i rfrf 5^ 3_ / 8 A / 2 



D 7 — 15- d^ 5 * D 3 ^ 5 ' D 5 ' 

~ ~ 35 * d^ 2 7 ' D 5 "+" 7 - D 7 ' 



W 



/ 3 



r 



D U "~ "~ 63 * dvf/ 2 — 9 # D 7 + 9 - D 9 * 

By means of the two first of these last formulas, we can exterminate yj- fi and 

f 3 ( 2 ) ( 3 ) 

fr^ from the values of B and B ; which being done, we shall find 

XJq 

t* (2) _(_L _L _L dd Vo ) >z + h f_ 



35 



1 \ / f 

r (3) _ f Izi J_ 9 + qk __5o I , ( 4&' + 9* + 9 / J_ ^ ^D7 



24 * D 12 * d^ 2 / ^ \ 24 ' D 3 12 * d$*' 

And, in like manner, by means of all the four formulas, we may exterminate 

f% fS f4 f5 (4) (5} 

IT* 5 IF"* IT 5 anc * D 11 ' ^ rom ^ e ex P ress * ons of B and B ; and thus are 
obtained these values, 

R ( 4 )_ 3 A 8 + 8 & ~ 7 1 , SF + 8F+8H6 / 
13 "~ "" 96' • D + 24 * D 3 



1 /" /" 



2 



5 



12^ + 32^ + 77 D J Dp 3 2. D o 

96 • tf^ 2 — 6 * rf*« ~32 ' dV ' 

MDCCCXXXIII. 4 E 
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(5) _ 24 A 3 4- 75h* 4- mh~ 55 J^ 96 h* + 300 F 4-341 If + 200 A + 159 f^ 
B — ~~ 960 * D Q + 960 * D^ 3 

1 f 



96 7z 3 + 300 A 2 + 368 h 4- 725 D 3 12 A a +££1 D 3 

960 * dty 2 960 * drj; 



2 



f2 f3 

25-4h a D 5 1 aa D 



7 



320 * d\[/~ 3 64 * d$* ' 

If we take the second differentials relatively to if/ in the first two of the same 

f 2 fs 

four formulas, we shall obtain these values of ■ d « 2 Q - and • d . a ° , viz. 
/* i l f 

(I rl J J*-~- dd •■-—- dd-— 

D 5 _ 4 ^ D 1 flfl D 4A ^Do% 



d^/ 2 — 3 • tfvf/ 4 3 • d*« ^ 3 • d^« 

,7,/«/L_ ,/* ,v ^* JL fl fi — /J rl * 

aa B " _ 4 ^ D 3 32 A D Q 8A D 32& g -9 D 3 



d^ 3 ~~ 15 * d^ A 15 ' d<J>* 15 ' dV ^ 15 * dty* * 

which being substituted, the expressions of B and B will be as follow : 

dd~ , d 4 l 



•R (4) — 3 A 2 -f 8 A - 7 _1_ 3ft 9 + 8A 4 19 D 2 Dp 

B ~~ ~~ 96 ' D 24 • d$* +24* d^ 4 

3h s 4-SF4- 8A 4-6 / H4 dd IV 



B 



24 * D 3 24 ■ • d ^ * 

1 
(5) <J 24 A 3 4- 75 g 4- 92A-55 J^ _ 24 A 3 + 75A 2 +9U— 175 ^fi" 

960 * D n 240 * d^ 

1 



d 



4 




25 4- 4ft D 

^ 240 ' dv^ 

96 ft 4 4- 300 ft 3 + 34JJ1+ ° 00 h + 1 5 ^ / _ 7 h* + 25 k + 53 d d D Q 3 
+ \ 960 * D 3 240 * d V ' 




240 ' dty 



Thus the six quantities B , B \ &c, are all expressed in terms of yr- 
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and jjs, and the second and fourth differentials of the same two quantities, 

the coefficients being of easy computation, and depending" solely on the mean 
distances of the two planets. 

5. The problem to be solved is, to expand jy in a series of terms arranged 

according to the cosines of the arc %|/, or <p — <p' + a — </, and its successive 

l f 

multiples. Now as the quantities jy and jp and their differentials, are all 

susceptible of a like development, what is required will be accomplished 
merely by substituting in B , B , &c., the parts of jp jy$, and their differ- 
entials, multiplied by the cosine of the same multiple of %//, instead of the 
quantities themselves. The coefficients of the cosines in the developments in 

ct ct 

question are functions of —,- or — ; and, for the sake of convergency, we must 

choose that one of the two fractions which is less than unit. Supposing that 
a! is greater than a, and that n represents any odd number, we have, 



n — I 
2 



/ 






X 



Wherefore, assuming as usual, 



= iC + C v 'cos^ + C ^cos2^... + C ^cos*^..., 



2 



the part of jj-, and of its second and fourth differentials, which are multiplied 
by cos k 4> 9 will be respectively, 

cos k 4 / c (*) cos A j; ^ _ ^ 2 c (^ cos/rtfr x ^,4 c (/ ° . 

and the like parts of jp and its differentials, will be, 

a 2 \„°L l _ ?1 

— r~ X C 3 cos&^, —zr—X -k 2 C 3 cos & ^, — 7— X A 4 C 3 cos&^. 

4 e 2 
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It remains then to substitute these values in B , B , &c, in place of jy-, jr-^ 

and their differentials : but, in order to obtain the formulas most convenient 
for calculation., it will be requisite to express the coefficient of the develop- 

ment of yts by means of those of the development of f=r , 
If we write a for -r, and put 

V = sj 1 +a 2 -2« cos ^, 
we have this formula, 

d - 

V a sin v[/ 



3 ? 



d4f — V* 

1 1 

and if we substitute the expansions of y and y 3> and equate the coefficients of 

sin k $ 9 the result will be • 

(*>■ (*-0 ~ (ft- 1) 

2 yt C x =:. # C 3 — # C 3 

We have likewise this identical expression, 

1 V 2 1 + O? %CL COS vj/ 



V — ys — v 3 "" Y 3, 

1 1 

and by substituting the expansions of y and y*, and equating the coefficients 
of cos ft -^ we g* et 

Ci = (1 + # 2 ) C 3 — a (J 3 — « C 3 

By combining the two formulas, we deduce 

2 a C 3 ( * + 1} .= (1 + a 2 ) C 3 W - (2 A; + I) C^ 

2 a C 3 *~ * = (1 + a 2 ) C 3 W - (2 & - 1) C/* 1 . 
Change & for & — 1 in the first of these formulas, and for A; + 1 in the second, 
and the two following formulas will be obtained, 

2 a C 3 W = (1 + <* 2 ) C 3 *~ ' - (2 A - 1) c/*" 

2*C 3 =(1 + « 2 )C 3 + (2&+ 1)C X . 

By combining the four last formulas in different ways, different values of C 3 
will be founds which may be used at pleasure, viz. 
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(i-rf»e i »-p* + i).{]i£q w -£ ? .c I * + ' ) } 

(1 -«*)C 3 «> = <»*_,) .{-^ C,"-"- 1±5X3,«} 

(l-^)C3 m = i^. T iV,.(c 1 1 - 1 -C +, ). 

In order to shorten expressions, and to bring the formulas to a form convenient 

for calculation, let 

u = tan | ; 
then, 

1 + a 9 , 1 2u sin , . , 

1""~ ' 2 — • #£> •— • TTT^fli "i 9 •"— "" /\ •-»— • it bill P • 

— or cos 0' 1 — a 4 cos 

further, assume the new symbol a , the value of which may be calculated by 
any one of these formulas, viz. 

W = (2A+l).(c i i -sm^C 1 (i + 1) ), 



a 



a 






then, 



,Bihtf(C 1 *"' 1 -C 1 * + 1 ); 



{\-~c£)C^ =zh.^\ 



f 
Using this value, the parts of jo and its differentials, multiplied by cos h $ 

will be 

^X k«^ x -K»*«\ ZZp X Aft»a*. 

These, as well as the like parts of D0 and its differentials, must be substituted 
in B , B , &c. As there is no longer occasion to refer to any development 
but one, C may be written for C x . In substituting we omit the common 

cos Jc ll/ 

factor — -, — , and put 
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for the expressions into which 



are changed, so that, 



-R (0) -R (1) -R (2) ™ (5) 



b ik) = c w 



1 — a ""■" a ' 



»? = (^ + 9 -4#-* • » ) ■ c w + r + r +iM + » • .», 

« / 3/z 2 -j~ 8^~7 , 3& g + 8&-M9 h2 , #A ^M 
, /3 A 4 + 8 h 3 + 8 A 3 + 6 h . A 3 + 4 A 79 \ (&) 

+ V — £4 — + -^ r ~.& 2 ).a •, 

(*) _ / 24 A 3 + 75 & g + mh -55 24F+ 75 k* ±9±h± 175 , 2 
»5 — V" - . 960 + £40 * * 

25 + 4k -A ^(*) 



i gj+1* jlA 

~ 240 ' n / * 



c 



, (96& + 300W + 34lh s + 200h2+\59h 7h 3 + Q5h* + 53h f9 , ££ 4 ^ (*) 

And, these symbols being established, the part of jy we are in search of, that 
is, the part multiplied by cos h Qp — <p ! + <r — <r% will be thus expressed, 

cos k (<p — $ -f <r — <r' 
f (A) (*) (*) (*). (*) (ft) 

) h + ^i * * l + h • * 2 + ^3 • ^ 3 + *4 * s' 4 + h • ^ 5 

* (A) (*) (*) (^) 

<-A6 .* — A&j .2^— A£ 2 .3s' 2 s— A6 3 .4A-AJ 4 . 5^ 4 ^ 

(ft) (^) (*) (ft) 

+ A 2 V • ^ + A 2 &x ,3^^+A 2 & 2 .6^^+ A 2 6 3 . 10*V 

(ft) (*) (ft) 

- A 3 b Q . s 3 - A 3 Ji .4^- A 3 6 2 .10^3 

(ft) (ft) 

+ A 4 ft -* 4 +A 4 &i .5^i^ 4 

(ft) 

""~" ZA c/a • t3 • 
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By the foregoing analysis the problem is reduced to a series of terms com- 
prehended in the formula, 

i ( Ic) v * 

± A h. t X ft / s l cos k (<p — £>' + <r — d) 

i f and i representing all positive numbers, zero included, of which the sum i f -f i 

does not exceed 5 ; j3 standing for the coefficient of / s* in (,?' + s) 1 % ; and the 
upper or lower sign taking place according as the index of A is even or odd. 
Although all the coefficients are supposed to be deduced from six of them, for 
which alone direct expressions are given, it is nevertheless obvious that an in- 
dependent formula may be found by which any proposed coefficient may be 

separately computed : namely, by substituting the values of b , b x , &c, in 
the expressions of the several finite differences. 

Although there can be no difficulty in applying the formulas, it may not be 
improper, for the sake of illustration, to take an example. We shall choose 
the instance of 

Venus and the Earth. 

In this case d will represent the mean distance of the Earth, and a that of 
Venus : wherefore, 

-^ = cc = 07323332 = tan \&, 0^71° 45' 32", 
h = — X —, = 3-1947202, log. sin 6 = 9*9776083. 

The method likewise requires that there be previously calculated a sufficient 
number of the first coefficients in the expansion, 

1 (0)0) (2) (k) 

^T + a ^c 2acos4/ = l c +C cos^+C cob2*....Ccos**...., 

which is accomplished by rules amply detailed in all the treatises on Physical 
Astronomy. We next deduce the value of a , viz. 

a = (2 k + 1) . \C — sm 6 C ) : 
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and, with the numbers thus founds we have generally for any proposed value 

of k, 

(k) (k) 



b x ik) = %C W + 1 -59736. a W 

b 2 =V-4+TJ' C +4-148919. a 



<*) 



b 3 ° £) — {- 0-008113 + 1-282453 . A 2 } . C W 
+ {10*459698 + 0-266227 k 2 } . a W 

fe 4 = j ~ 0-512255 + 3-132352 * 2 + |- } . C 

(ft) 
+ {28-09035 + 0-957713 A 2 } • a 

(ft) (*) 

\ ={- 1-861383 + 8-390554^+ 0-157412 A^}.C 

(k) 
+ {80-06813 + 2-719660 k 2 + -0133113 &} . a 

Nothing ean be more easy than the computation of these six quantities for any 
assigned value of k. The other quantities, being the several orders of the finite 
differences of the first six, will be known when these have been computed: but 
it will be convenient and will often save much calculation to have an inde- 
pendent formula for every coefficient separately. 

- A b Q ik) = £ C W _- 1-59736 . a W 

A , (ft) / 1 &\ ~(*) (k) 

- Ab x — (^ — -$) C — 2-551559. a w 

(ft) (ft) 

- A b 2 = {0-2581 13 - 0-782453 A 2 } . C 

(ft) 
— {6-31078 + 0-266227 k 2 } . a 



(ft) f mi (k) 

> 3 = j 0-504142 — 1-849899 i 2 -i[C 



{17*63065 + 0-691486^} . a W 
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(ft) (ft) 

Ab t = {1-34913 — 5-258202 A 2 — 0' 11 5745 A 4 } C 



— {51-97778 + 1-761947 k 2 + -01331 13 &*}.a (/ ° 
AH « = (1 + 1) C W + 0-954199. a W 

A 2 ^ ( '°= { — -0081 13 + -282453 A 2 } . C 
+ {3*75922 + -266227 k 2 } . a W 

A 2 & 2 (/<:) = { — -24603 + 1-067446 A 2 + -041667 &} . C 

(ft) 
+ {11-31987 + 0-425259 k 2 } . a 

(ft) (ft) 

A 2 b 3 — { - -84499 + 3*408303 A 2 + -074078 A*} . C 

(*) 
+ {34-34714 + 1*070461 A 2 + '01331 13 £*} a 

— A 3 6 W = {0-25811 + -217547 A 2 } . C (k) 

— {2-80502 + -266227 A: 2 } . a W 

— A3 & X W = {-23792 — '784993 A 2 — -041667 A 4 } . C W 

— {7-56065 + -159032 A 2 } . a W 

— A3& 2 W = {-59896 — 2-340857 k 2 — -03241 1 A 4 } . C W 

— {23-02727 + -645158 Jc 2 + -013113 A 4 } . a W 

(ft) (ft) 

A 4 V = {'020197 + 1-002540 k 2 + '041667 A 4 } C 

+ {4-75563 — -107195 A; 2 } . a W 

Ik) (k) 

A 4 ^ = {— -36104 + 1-555864 k 2 — -009256 A 4 } . C 
+ {15-46662 + -486126 A; 2 + -0133113 A 4 } .a 

(k) (k) 

— A 5 V = {"38124 — -553324 k 2 + '050923 A 4 } C 

(k) 

— {1071099 + -593321 k 2 + '0133113 ¥} . a . 

According to these formulas the calculation of any of the coefficients is an 
easy arithmetical process. In order to be able to compute, in any term of the 
mdcccxxxiii. 4 F 



582 MR. IVORY ON THE DEVELOPMENT 

expansion of ^ , the part of any proposed order with respect to e and e', no- 
thing more is wanting than a method for reducing 

/3 / S COS (<p — $ + (T — d) 

to a series of simple cosines. On this point some observations will be offered 
below : but, without proceeding further, there is no difficulty in this respect, 
when we confine our views to quantities not passing the second or third order, 
which comprehends all the perturbations useful in astronomy, except the 
inequalities of very long periods. 

If we combine the planets two and two, as is done in the sixth chapter of 
the sixth book of the M6canique Celeste, and for every two planets construct a 
formula such as is exhibited above for Venus and the Earth, limiting the extent 
of the calculation according to the nature of the case, the theory of the pla- 
netary disturbances would be rendered more accessible, and would be freed 
from the tedious and disgusting labour which has rendered astronomers averse 
from cultivating this branch of their science. 

6. The expression, 

s 1 s* cos k (<p — <p' + a — </) 

may be reduced to this form, 

{i i' i i' 1 

(& cos k a) (V cos k J) + (s sin k <r) (/ sin k <r') ? 

(s sin k(s) ($' cos k<r f ) — (s cos &<r) (s' sin&<r') k 

Here the quantities within the brackets are serieses, of which each contains 
the mean anomaly of only one planet, and comes under one or other of eleven 
different forms, namely, the six, 

cos k <r, s cos k <r, s 2 cos k <r, ...... $* cos k <r, 

and the five, 

sin k cr, s sin k or, s 2 sin k <r, . . . . . . s A sin k <r, 

the quantity $* sin k <r, which is of the sixth order, being omitted. None of these 
serieses, when carried to quantities of the fifth order inclusively, and arranged 
according to the sines and cosines of the multiples of the mean anomaly, con- 
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tain more than six terms, and they thus afford a ready way of calculating the 
part of 

i' i 

$' s cos k (jp — (p 1 + * — d) 



of any proposed order with respect to e and e f . 
These serieses are as follow : 

Cos&* = (l - k 2 e 2 - fo& e* + ^f) 



f L Jfi e 2 + L k 2 e 4 _ ^ & e ^ , e C os im 



+ fJLtf _ ?* ^4 e 2 _ £ &4 e 2^ e 3 cos 3 pi 



+ (ls +T92^j.e 4 cos4^ 
. /5A* , 7F\ - 



£* 9 . 5Ae* ,„ „ . 13Fe* . k b e*> 



Sm Ar <r = ( 2 A: j- + Q(i — A 3 e 2 + -gg- + -g~ j • e sm (j 

+ (T - ~£3 -4-) .e2sin2ffc 

, /ISA 43 7 9 F 307F.€» J^A , . „ 

. /1097 7 . 179,- . **\ 5 . . 

s cos A: <r = — -g- + -g~ 

/ Sf 8 . 5e 4 *V , 37 A 8 * 4 ; **A 

+ (l--F+T92 ^ + -96- + IT) • c cas f 6 

+ (-5 « ^ 2 e 2 j . e 2 cos 2 ft 

4 f2 
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3 45e 2 . k* 183. P e 2 F^ 



, /3 45 e* t &* 183. A* ^ A7***\ * 
+ (F-l28"+^ IsT^ ~ T") • «* C0S 3 ** 

(1 7 A a \ 
T + IT ) ' e4 cos 4 ^ 



125 , 475 TO , k 4 \ 



+ (384+3874^+24)- eCOs5 ^ 

$ sin # <r = — k I -j e 2 — 3§ ^ 4 £4™ ) . e sin f6 



(23 & £ \ 

, , /9 115. e 2 35Fe 2 \ , . „ 
+ A (¥ U6 - ~ ~48~ ) ' * Sin 3 ^ 

+ ( q" A + 4 8 &) . e 5 sin 5 ^ 
s 2 cos Ic cr = Y ~ — 

/e 2 e 4 , PA 

+ (a ) • e2 cos 2 p 

' /l. 9<? 2 7FA , _ 

+ (¥-Tti--T6-)- e3cos3 ^ 



- > 



+ (t + x) ' «* cos 4 



+ (i8 + 16 )'* C0S 5 ** 

. . , (Ice* , jfc«* ¥e\ 
■s 2 sin A ff = f -gj- 4-^g jg j . e sin f& 

g— . e 2 sin 2 [a 

+ ( T - "6 la") • ^ sm 3 ** 
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+ (l¥" + 12) • e Sm 5 <* 



3 



a 3 cos k & = — -g- e 4 



(8 8 <? 4 \ 

-j e 2 — -jq + & 2 e 4 J . e cos \h 



+ (t~ 32 + -§") • c 3 cos 3 ^ 



3 



-f -g . e 4 cos 4 ^ 
+ /q3 + -g-j . e 5 cos 5 ^ 
s 3 sm & <r = — -jg- . e sin ju + -£- . e 2 sin 2 /£ + -^- . e d sin 3 ^ 

h \1 h 

+ — . e 4 sin 4 ^ + -g^- . e 5 sin 5 ^ 

3 g 4 £ 2 £ 2 

£ 4 cos&<r = g-e? 4 — -5 .e?eos^ + -g«€? 2 cos 2 1&+--£. .e3 COS 3^ 

1 1 

4- -g . e 4 cos 4 ^ + "4 • ^ cos 5 ^ 

o sm^= — . e sm ^ H — . c 3 sm 3 ^ + — . e 5 sm 5 /*, 

5 e 4 5e 2 1 

♦s 5 cos k <r = -g- ,eco${jb + -jq- . e 3 cos 3 f* + ~yg e 5 cos 5 ^ 
The use of these serieses is obvious. If we wish to determine all the argu- 

a? x f 

ments of the order e e in the term multiplied by cos k (<p — qf + # — </), of 

1 
the development of ^ ^ and of representing two positive integer numbers,, zero 

included^ of which the sum does not exceed 5 ; we first set aside all the parts of 
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the expression in which i is greater than x, or i greater than x\ for these con- 

x af 

tain no quantities of the order e e' . Let 

i Qc) i V 

A b X |3 s s cos k (<p — <p' + <* — *') 
be one of the remaining parts, in which i is not greater than x, nor i' greater 



than x; in the serieses s cos k a and s sin k <r, take the terms A e cos x p and 

ir i f i x* 

Be sin x p; and, in s' cos k a and s' sin k </, take the terms A' e cos a?' y! 

x' x x' 

and B / e' sin # f//; then all the arguments of the order e e will be contained 
in the formula 

i (k) x x* 

A b., X |3 e e' X 

cos k(<p — <p') . {A A' cos x (ii cos x /// + B B' sin x p sin a?' p'} 
— sin& (<p — p') . {A'B sin# ^ cos^'f-o' — AB'cos<r ^sino?'^'}. 

And thus are computed all the quantities of the order e e' x in any term of the 
development of jr. The procedure is exactly the same with that followed in 
the expansion of the first part of the disturbing function in § 2. All the argu- 
ments of the order / e' are only four, comprehended in the formula 

jfc (<p — <p') + x f/> + x / (jJ, 

which represents all the different ways of combining x \h and x f p, by addition 
and subtraction, with k (<p — <p'). Hence it is easy to ascertain the orders of 

the same argument, or the relative magnitude of its coefficients, as it recurs 

i 
in different terms of the development of jj. 

7. We have next to consider the parts of the disturbing function which de- 
pend on the inclination of the orbits. 

1 T T* 

Development of sin 2 -~r I . ^ . q. 
Omitting quantities above the third order with respect to s and / in the ex- 
pression of jy found in § 3, we have 



u 
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J (0) C 1 ) ( 2 ) (3) 

^=B +B V+B V 2 +B V* 

(0) (1) (2) 

-AB *-AB 2n'-AB 3s*' 2 

(°) (0 

+ A 2 B V + A 2 B 3 s 2 s' 

(o) 

~A 3 B A 

(o) (i) 1 

As the values of B , B , &c, contains- and its powers, they may be re- 
garded as functions of cos %£: taking the differentials with regard to this vari- 
able, we get, 

dB ^ a a! a a 1 f a a! tv 

d"(cos*) — / I s • iy + 2 • W )~ f ' 

rfB^ 5 _ c t°t( < ^JL l Azljj Z 2 _l J£ /i i_ ^ tv (2) 
rf(cos"4i) ~ / y 8 b a + 4 • D 5 + 8 D 7 y— / x B * 

(3) / 

^ _ aaf__ I _5 /_ £7-18^ /* 45 — 30 h f*_ 35 f* 
d (cos $) — / X y 16 " B 3 + 16 ' b 5 + 16 * D 7 + T6 ' D 

= yXB' . 

i 
If we now differentiate the above formula for g, and observe that 

a a ' ^ a ' 1 / — tan ^ 

~J~ — V a ~ a 2 = -Q+jh 2 — 1 — — g-, 

the result, when the new symbols are introduced, will be as follows : 

rr ; tan. Q 



9 



D 3 _ 2 



f (°) (1) (2) (3) 

X | B A + B A V + B' V 2 +B A V 



(0) (1) (2) 

AB' *-AB' ,2^-AB' 8*^ 

,(o) (1) 

+ A 2 B' V + A 2 B' 3 a 2 a 1 

(°) 
-A 3 B' V. 
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(O) (1) 

The values of B' , B' , &c, may be reduced to expressions that contain only 

i / 

p- and j^t and their second and fourth differentials, by the same process that 

(0) (1) 

was applied to B , B , &c, in § 4: and in this manner we obtain. 



(0) 

B ; = 


""Do 3 


(!) 


_ 1 1 


TV 


3i + 2 


Jl> 


*^fc 



rfi- * 



" ^ tf* 9 "T" 2"^ 'Do 3 ' 

• ^ - (»* + *) •-# 

4-^3 A 2 + 2 A - — ^ -^o- - — r5° 3 - 

i i 

d d -— r/ 4 — - - 

« _ 24 y + 27 A + 1 J_ 46 A 9 + 54 A + 1 D J_ Dp 

13 _ ~~ 24 • D 12 * d^ + 3 * d^ 4 

96A 8 +1 08A 9 -23A--33 / 10^ + 9 ^~WL 

The term multiplied by cos & -^ in the expansion of --gr, will now be ob 

tained by substituting in B' , B , &c, the parts of-jj-, -j-pr, and their differ- 
entials, multiplied by the same cosine. Let 

cos k \J/ , ,( 7 cos 7c 4/ ,/( /f ) cos#\}/ , , W cos&vl; _ ,(*) 

denote the expressions into which 

BAN / __ *\ / "V*\/^ ~W*% / 

are changed by the substitutions mentioned ; the quantity sought, or the part 

of the expansion of sin 2 |lx -gr X 9 multiplied by cos & ^, will be thus ex- 
pressed, 
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sin 2 \ I . cos Jc (<p — <p' + <r — or 1 ) x q 



a' X 



tan 9 f .,(*) , .,(*) , , „(*) ,, , I( (*) 



z> 



b + b\ $ + &' 2 s' % + b' s *' 3 

(*) (70 (*) 

— A6' s — Ab\ 2ssf—A6\ 3ss 2 

(Ic) He) 

+ Ai', s 2 + A 2 //' 3 s 2 / 

and the values of b' Q , b\ , &e., will be as follow: 

/; = h ,a , 

b \ = (- i + 2 A 2 ) C + — 2 • a 

tf, ==(3A2).-^— i.C + (3A 3 + 2A 2 f A+^.A 2 ).a 

,,(*) f 24 A 2 4- 27 A 4- 1 . 48 A 2 + 54 A +1 7a . AM ~(*) 

6 3 =| ^™ + la •* +"3~j- c 

. f96A 4 + 108 A 3 -23 A 2 - 33 A 10A 2 +9A 1 (*) 
+ ) 24 + 12" A /- a ' 

By taking the several orders of the finite differences of these four quantities, 
the other six coefficients may be numerically computed, or an independent 
formula may be found for calculating each of the six separately. The expan- 
sion of this part of the disturbing function is therefore reduced to the ex- 
pressing of 

if i 

j8 s' $ COS A (<p — • <p + or — a) X q 

in a series of simple cosines, 1 and i being any positive numbers, zero included, 

if i 

the sum of which, 1 + h does not exceed 3, and (3 being" the coefficient of / s 
in ($' + s) . By substituting what q stands for, we have, 



1 1 



s 1 s cos k (<p — <p f -{-ry — (r f ) X q 

i f i 
\ d s cos ( A — 1) (<p — <p' + fi — • </) 



X I 



+ \ d s cos (A + 1) (<? — <?' + * — 
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¥ I 



- I s s cos (ft — 1 . — <p ft + 1 * p* ~ 2 v + ft - 1 . (T — ft + 1 . *') 






— fc^* cos(ft + l.<p — ft- l.p'— 2* + ft + 1.* — A— l.<r r ). 

As the arguments in the expansions of these four terms are different, each 
term must be computed separately, for which purpose the method in § 6 may 
be used. Thus 

/s cos(ft~^l.p— ft+ 1 . p' — 2* + F^~l. <r- ft+1 .<r r ) 
■= cos (ft^"l . <p - ft+l.<p'~2?). | (/cos (ft— 1).*) (/eos(ft-f 1)/) 

+ (/sin (ft - 1) <r).(/sin (ft + 1) *) I 
~~ sin ( F^l . <p — ft+1 . ^ r - 2 • { (/ sin (ft — 1) <r ) . (/ cos (ft + 1) tf 

- (/ COS (A - 1) *) . (/ sin (ft + 1) </) | ; 

and the quantities within the brackets being known serieses, the part of the 
expansion of any order / e 1 , is readily obtained. 
8. One part of the disturbing function yet remains to be considered. 

Expandon of sin 4 £ I X f . -jjr- X q\ 

In the expression of -jjr found in § 7> leave out quantities of the second and 
third order with respect to s and $\ and we shall have, 

' rr 1 tan. 8 f 1 w (p) (i) 

i-D 5 = ~Y" X |£B' +IB 1 .e'cos^' 



Now let 



o 
-| AB' . ecos^. 

i rf g(0) — i^ *-ZL _ ^ B " (0) 

a # d(cos^) ~ f m * D Q 5 — f X 

^ ^(costf/) - / V*-D 5 + 4 * D V~ / X ^ ; 
and differentiate the formula for | . ^ relatively to the variable cos 4s intro- 
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ducing the new symbols, and observing that -y- = —^- , the result will be, 



r 2 r m tan 2 

3 



^x{b'' ( %B'' (1 Vcos^' 







- AB ,; e cos p. 

(0) (1) 1 f 

The quantities B ' and B" , when expressed as before in terms of g- and -d^ 



and their differentials, will be as follow : 
B -|. Do 5- - 3f Do - 2 -^2 + J *-D % 



d*< 



G — 4 ' D 5 "*" 4 # D 7 



LA /• . j j J 



- (8 A +1). I ±."5- + *^ J- (8A 2 + A-|).^3 

\ ° ~dV / ° 

Proceeding as before, let 

COS k \{f _ ,,(&) , cos & \f/ 7 .. W 

— ^r 1 X b\ , and — ^r~ X b\ , 
represent the expressions into which 



d d 



D3 
n 



ay* 



B /,(0) and B ,,(1) 



1 / 

are changed, when we substitute, instead of g™ and-jjs and their differentials, 

the parts of these quantities multiplied by cos k $ : then the part multiplied 

by cos k ^ in the expansion of sin 4 | I X § -jjr- X <f, will be thus expressed: 



>in 4 \ I X cos & (<J> — <p' + cr --or') x 5 



sin* -1 1 x cos/t ($ — p' + cr —or'; x £ 2 

x 



4 



X <j 6" + 6 \ . <? cos yd 



Ab n e cos ^; 



and the values of b n and b x will be as follow : 






4 G 2 
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The coefficients for any proposed value of k may now be readily computed ; 
and nothing is wanting in this part of the problem but to reduce the factor^ 

cos A (p — p' + «■ — ff ). < cos (<p — • p + * — *0 — ■ cos (<p + p + <y + $ — 2 ') r , 

to the form, 

M + M'xe sin /* + M" X e' sin /*/, 

and then to multiply by it, retaining in the product only quantities of the first 
order relative to e and e. Now in this there is no difficulty : but the complete 
development of the expression would be bulky, and would contain many argu- 
ments, the greater part of which are insignificant and useless, although some 
of them deserve attention in particular researches. Leaving astronomers to 
select from the general expression the arguments which may suit their pur- 
pose, we shall here close what we had to offer respecting the development of 
the disturbing function, without adding to the length of this paper by any 
application of what we have written. 
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